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Summary 

This  investigation  aims  at  the  elastostatic  field  near  the  edges 
(tips)  of  a plane  crack  of  finite  width  in  an  all-around  infinite  body,  which 
at  infinity  — is  subjected  to  a state  of  simple  shear  parallel  to  the  crack 
faces.  The  analysis  is  carried  out  within  the  fully  nonlinear  equilibrium 
theory  of  homogeneous  and  isotropic,  incompressible  elastic  solids. 
Further,  the  particular  constitutive  law  employed  here  gives  rise  to  a 
loss  of  ellipticity  of  the  governing  displacement  equation  of  equilibrium 
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in  the  presence  of  sufficiently  severe  anti-plane  shear  deformations. 

The  study  reported  in  this  paper  is  asymptotic  in  the  sense  that 
the  actual  crack  is  replaced  by  a semi-infinite  one,  while  the  far  field  is 
required  to  match  the  elastostatic  field  predicted  near  the  crack  tips  by 
the  linearized  theory  for  a crack  of  finite  width.  The  ensuing  global 
boundary- value  problem  thus  characterizes  the  local  state  of  affairs  in 
the  vicinity  of  a crack-tip,  provided  the  amount  of  shear  applied  at 
infinity  is  suitably  small. 

The  results  communicated  in  this  paper  were  obtained  in  the  course  of 
an  investigation  supported  by  Contract  N00014-75-C-0196  with  the  Office 
of  Naval  Research  in  Washington,  D.  C. 
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An  explicit  exact  solution  to  this  problem,  which  is  deduced  with 
the  aid  of  the  hodograph  method,  exhibits  finite  shear  stresses  at  the  tips 
of  the  crack,  but  involves  two  symmetrically  located  lines  of  displacement- 
gradient  and  stress  discontinuity  issuing  from  each  crack-tip. 

Introduction 

In  a number  of  investigations  conducted  during  the  past  few  years 
we  have  explored  two  distinct  types  of  inherently  nonlinear  singular 
problems  in  finite  elastostatics.  The  first  of  these  categories  pertains 
to  geometrically  induced  point  or  line  singularities,  such  as  those  arising 
in  connection  with  crack  and  punch  problems,  for  which  the  linearized 
theory  of  ealsticity  — in  conflict  with  its  underlying  approximative  as- 
sumptions — predicts  locally  unbounded  deformation  gradients.  The  second 
category,  which  has  no  counterpart  in  linear  elastostatics,  concerns  the 
emergence  of  lines  or  surfaces  across  which  the  deformation  gradients 
suffer  finite  jump  discontinuities.  This  kind  of  singular  behavior  becomes 
possible  only  if  the  elastic  material  is  such  as  to  admit  a failure  of  el- 
lipticity  in  the  corresponding  displacement  equations  of  equilibrium  at 
sufficiently  severe  deformations. 

The  results  obtained  in  these  investigations  are  summarized  in  a 
comprehensive  survey  [l],  which  also  contains  references  to  the  related 
background  literature*.  The  present  paper  extends,  and  combines  strands 
from,  both  strings  of  studies  alluded  to  above.  It  deals  with  the  elasto- 
static  field  near  the  tip  of  a crack  under  conditions  of  anti-plane  shear  in 
a constitutive  setting  within  the  fully  nonlinear  equilibrium  theory  that 


*In  this  paper  we  shall  cite  explicitly  only  previous  work  that  is  of  im- 
mediate relevance  to  our  current  purpose. 
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leads  to  a loss  of  ellipticity  and  the  concomitant  appearance  of  lines  of 
deformation-gradient  discontinuity  ("elastostatic  shocks")  in  the  vicinity 
of  the  crack-tips. 

The  particular  problem  treated  here  is  suggested  by  and  complements 
an  asymptotic  study  reported  in  [2],  We  consider  an  all-around  infinite  — 
homogeneous  and  isotropic  — incompressible  elastic  body  with  a plane 
crack  of  constant  width  and  infinite  length.  At  infinity  the  body  is  sub- 
jected to  a state  of  simple  shear,  parallel  to  the  faces  of  the  crack.  In 
[2]  the  deformations  and  stresses  ensuing  near  the  tips  (edges)  of  the 
crack  were  analyzed  for  a class  of  such  materials  whose  elastic  potential 
(strain-energy  density)  is  completely  determined  by  the  shear-stress 
response  to  a simple- shear  deformation.  Moreover,  the  analysis  carried 
out  in  [2]  was  confined  to  circumstances  in  which  the  induced  shear  stress 
is  a steadily  increasing  function  of  the  amount  of  simple  shear  applied. 

In  this  instance  the  displacement  equations  of  equilibrium  are  found  to  be 
elliptic  at  all  field  points  and  at  every  solution  corresponding  to  anti-plane 
shear.  In  contrast,  ellipticity  is  lost  at  severe  enough  anti-plane  shear 
deformations  if  the  response  curve  in  simple  shear  undergoes  at  least 
one  reversal  in  slope. 

The  case  in  which  the  response  curve  in  simple  shear  has  a 
positive  slope  throughout  but  rises  monotonically  to  a finite  shear  stress  — 
as  the  amount  of  shear  grows  beyond  bounds  — thus  signals  an  impending 
breakdown  of  ellipticity,  which  would  occur  if  the  stress  were  to  increase 
steadily  to  a maximum  value  at  a finite  amount  of  shear  and  thereafter 
declined  monotonically.  This  troublesome  transition  case  is  still  fully 
discussed  in  [2]  with  the  aid  of  the  hodograph  method  and  on  the  basis  of 
an  asymptotic  scheme  in  which  the  crack  of  finite  width  is  replaced  by  one 
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of  semi-infinite  width,  while  the  far  field  is  required  to  match  the  elasto- 
static  field  near  the  crack-tip  predicted  by  the  solution  to  the  original 
problem  according  to  the  linear  theory. 

The  results  deduced  in  [2]  for  the  limiting  elliptic  case  described 
above,  though  of  unlimited  smoothness  in  the  interior  of  the  body,  display 
certain  features  that  suggest  the  emergence  of  discontinuities  in  the 
deformation  gradients  once  ellipticity  is  lost.  In  order  to  explore  this 
eventuality  by  means  of  a pilot  example  we  presuppose  here  a specific 
response  curve  in  simple  shear  that  rises  linearly  to  a maximum  shear 
stress  and  then  declines  steadily  as  the  amount  of  shear  increases.  The 
special  constitutive  assumption  thus  introduced  enables  us  to  establish 
on  the  basis  of  the  hodograph  method  an  explicit  exact  solution  to  the 
crack  problem  under  present  consideration.  The  solution  thus  obtained 
involves  two  lines  of  displacement- gradient  and  stress  discontinuity  that 
are  symmetrically  located  with  respect  to  the  axis  of  the  crack,  issue 
from  the  crack-tip,  and  terminate  in  the  interior  of  the  body. 

1 . Preliminaries.  Finite  anti-plane  shear  for  incompressible  elastic 
solids.  Formulation  of  crack  problem. 

In  this  section  we  cite  from  [2]  certain  prerequisites  needed  here. 
Suppose  the  interior  of  a body  occupies  an  open  region  ft  in  an  undeformed 
configuration  and  consider  a deformation 

y =y(x)  =x  + u(x)  for  all  x€ft  (1.1)* 

that  maps  ft  onto  its  deformation  image  ft  . Here  x is  the  position 
vector  of  a material  point  in  ft,  ^(x)  its  position  vector  in  ft  , and  ^u(x) 



Boldface  letters  are  used  to  denote  vectors  or  tensors. 


r 


its  displacement  vector.  We  assume  the  mapping  (1.  1)  to  be  suitably 


smooth  and  invertible.  Let 


F=Vy=l+Vu  on  ft  , 


(1.2) 


so  that  F(x)  is  the  deformation-gradient  tensor  at  x.  For  an  incompres- 
sible  medium  the  deformation  must  be  locally  volume  preserving  and  thus 


J = det  F = 1 on  ft 


The  left  deformation  tensor  associated  with  (1.1)  is  given  by 


G = FF 


and  has  the  fundamental  scalar  invariants 


(1.3) 


(1.4)' 


I.  = trG  , I-=-ktrG)2-tr(G2)]  , I,=  detG  = J2=  1 . (1.5) 

sjc 

Next,  let  T be  the  actual  (Cauchy)  stress  tensor  field  on  ft  and 


0 the  corresponding  nominal  (Piola)  stress  field  on  ft,  whence 

T T -1 

T=CTFi  , 0 = T(FX) 


(1.6)' 


The  balance  of  linear  momentum,  in  the  absence  of  body  forces,  leads  to 


the  equilibrium  equations 


5)C 

divT=0  on  ft  or  div(J  = 0 on  ft  . 


(1.7) 


Suppose  further  the  incompressible  medium  at  hand  is  a homo- 
geneous, isotropic,  elastic  solid  that  possesses  a strain-energy  density 
W per  unit  undeformed  volume.  In  this  instance  W = Wfl^,^)  and  the 

^We  write  1 for  the  idem-tensor. 

2 

The  superscript  T indicates  transposition. 

3 — 1 

If  A is  a nonsingular  tensor,  A stands  for  its  inverse. 


■ ' 


appropriate  constitutive  law  takes  the  equivalent  alternative  forms: 


S=2[lfe  + |r<I1A-5)F]-P(FT)-1  . 


(1.8) 

(1.9) 


where  p is  the  arbitrary  hydrostatic  pressure  required  to  accommodate 
the  kinematical  constraint  (1.3).  The  stress-deformation  relation  (1.8) 
assures  the  symmetry  of  £,  but  £ is  in  general  not  symmetric. 

J{C  # 

If  S is  a surface  in  ft  and  S its  deformation  image  in  ft  , 

one  has 


Tn  = 0 on  S if  and  only  if  an  = 0 on  S , 


(1.  10) 


# 5JC 

where  n and  n'  are  the  respective  unit  normal  vectors  of  S and  S . 
Thus,  in  particular,  a portion  of  the  boundary  of  ft  is  free  of  actual 
tractions  provided  the  nominal  traction  vanishes  on  its  antecedent  within 

the  boundary  ft.  Further  the  actual  traction  t = Tn*  is  found  to  be 

* 

continuous  across  S if  and  only  if  the  nominal  traction  s =an  is 
continuous  across  S. 

With  a view  toward  the  crack  problem  that  constitutes  our  objective 
we  now  assume  that  ft  is  the  exterior  of  a plane  infinite  strip  of  width 
2c.  Further,  we  introduce  rectangular  cartesian  coordinates  (Xj,x2,x^) 
in  such  a way  that  this  strip,  which  represents  the  crack,  lies  in  the  plane 
x2=0,  the  x^-axis  being  parallel  to  the  edges  of  the  crack  and  the  origin 
midway  between  the  edges.  Figure  1 shows  the  cross-section  & of  ft 
in  the  plane  x^  = 0.  Evidently,  the  plane  domain  f)  is  the  exterior  of  the 
line  segment  x2=  0,  -c^x^c  in  the  (x^,  x2)-plane.  We  confine  our 
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attention  now  to  the  case  in  which  the  body  is  subjected  at  infinity 

to  a simple  shear  parallel  to  the  crack  faces *.  The  displacement 

2 

components  are  thus  required  to  satisfy 

ua=o(l)  , u3=kx2+o(l)  as  x^x^-oo  (-oo<x3<oo)  , (1.11) 

in  which  the  positive  constant  k is  the  amount  of  the  applied  shear. 

Moreover,  the  deformed  faces  of  the  crack  are  taken  to  remain  free  of 
actual  tractions.  In  view  of  (1.  10),  this  assumption  yields  the  boundary 
conditions 

Oi2(xi,0±,x3)  = 0 (-cCXjCc,  -oo<x3<oo)  (1.12) 

for  the  components  a.,  of  the  nominal  stress  field.  If  ft  is  understood 
to  be  the  particular  domain  just  defined,  the  crack  problem  to  be  con- 
sidered accordingly  consists  in  finding  a deformation  (1.  1)  obeying  the 
incompressibility  condition  (1.3),  as  well  as  a hydrostatic  pressure  field, 
such  that  the  displacements  satisfy  (1.  11),  while  the  associated  nominal 
stress  field  (1.9)  satisfies  the  equilibrium  equation  (1.7)  and  meets  the 
boundary  conditions  (1.  12). 

At  this  stage  we  restrict  ourselves  to  the  class  of  incompressible 
elastic  solids  for  which  W is  independent  of  the  second  deformation 
invariant  and  set 

W = W(Ij)  for  all  Ij  s 3 , W(3)  = 0 . (1.  13)3 

^his  loading  case  is  known  as  "Mode  III"  in  the  terminology  of  fracture 

mechanics.  I i 

2 

Latin  and  Greek  subscripts  have  the  ranges  (1,2,3)  and  (1,  2), respectively. 

Summation  over  repeated  subscripts  is  taken  for  granted. 

3 

It  is  easily  seen  from  (1.4),  (1.  5)  that  1=3  at  the  undeformed  state  and 
that  otherwise  Ij>3.  1 


Such  a solid  has  a shear  modulus  at  infinitesimal  deformations  given  by 


* H = 2W'(3)  (1.14)1 

and  conforms  to  the  Baker- Ericksen  inequality,  provided 

W/(I1)>0  for  Il>3  . (1.15) 

In  particular,  for  a NeoHookean  material, 

W(Il)=2(Il"3)  • (1-16) 

For  materials  governed  by  the  strain-energy  density  (1.  13)  the 
foregoing  crack  problem  is  consistent  with  the  assumption  that  the  ensuing 
deformation  throughout  ft  is  one  of  anti-plane  shear  parallel  to  the  faces 
of  the  crack: 

u =0  , u =u(x.,x?)  on  ft  . (1.17) 

Indeed,  as  shown  in  [2],  (1.  17)  enables  one  to  reduce  the  original  three- 
dimensional  problem  to  the  following  subsidiary  two-dimensional  boundary 
value  problem  for  the  unknown  out-of-plane  displacement  u on  the  cross- 
sectional  domain  &: 

[W'(I.)u  J =0  on  £>  , (1.  18)2 

1 , Qr  9 Oc 

with 

Ij=3+  |Vu|2  , |7u|2=u>au>a  ; (1.  19) 

^Here  and  in  what  follows  W7(Ij)  stands  for  the  derivative  of  W with 
respect  to  its  argument  Ij. 

2 

Subscripts  preceded  by  a comma  indicate  partial  differentiation  with  respect 
• to  the  corresponding  material  cartesian  coordinate. 


i 


u(x, , x0)  = kx0+ o(  1)  as  x x —co  , (1.20) 

i c.  c.  a a 

u 2(x^»0±)  = 0 (-c<Xj<c)  . (1.21) 


The  hydrostatic  pressure  p = 2W/(I^)  on  f)  in  the  present  instance  and 
(1.17),  (1.  13)  together  with  (1. 8),  ( 1 . 9)  yield  the  stress-displacement 
relations 


T3a=Vr°3«='Vr2W'<Il>“.«  • 

(1.22) 

TaP=°ttP=a33=0  ' T33=  2w'<iiHvu| 

(1.23) 

with  1^  supplied  by  (1.19). 

We  now  specify  the  smoothness  demanded  of  the  solution  u to  the 
boundary- value  problem  characterized  by  (1.18),  (1.19),  (1.20),  (1.21). 

In  this  connection  we  suppose  that  W has  a continuous  first  derivative 
and  at  least  a piecewise  continuous  second  derivative  on  its  interval  of 
definition  [3,  co).  Further,  we  restrict  u to  be  continuous  and  to  possess 
piecewise  continuous  first  and  second  partial  derivatives  on  £;  in  addition, 
u is  to  be  bounded  on  every  finite  sub-domain  of  j6.  Thus,  in  particular, 
u must  remain  bounded  within  any  circle  of  finite  radius  centered  at 
either  crack-tip.  Next,  the  limits  Vu(Xj,0±)  are  to  exist  and  be  continuous 
for  -c<Xj<c.  Observe  that  the  preceding  smoothness  restrictions  allow 
for  the  possibility  of  finite  jump  discontinuities  in  Vu  — and  hence  also  in  the 
stresses  — across  curves  in  &.  If  J is  such  a curve,  however,  we 
stipulate  that  W^IjJSu/Sn  be  continuous  across  provided  n is  the 
unit  normal  vector  of  J and  9u/3n=Vu»n.  Let  s = On , so  that  s is 

r 

the  nominal  traction  vector  along  •/.  One  draws  from  (1.22),  (1.  23)  that 


-10- 


(1.24) 


The  foregoing  continuity  condition  therefore  assures  the  continuity  of 
the  Piola  tractions  across  the  cylindrical  surface  whose  trace  upon  f) 
is  J . Consequently  , the  Cauchy  tractions  are  continuous  across  the 
deformation  image  of  this  surface.  The  continuity  of  both  the  nominal 
and  the  actual  tractions,  in  turn,  is  a necessary  equilibrium  requirement. 

If  the  boundary  conditions  (1.21)  are  deleted  and  & is  momentarily 
taken  to  be  the  entire  (x^,  x^J-plane,  i.e.  in  the  absence  of  a crack,  the 
differential  equation  (1.  18),  subject  to  the  condition  (1.20)  at  infinity, 
clearly  admits  the  elementary  solution 


u(x1,x2)  = kx2  for  all  (x^x^eS  , (1.25) 

which  corresponds  to  an  undisturbed  homogeneous  field  of  simple  shear 
of  amount  k.  Also,  (1. 22),  (1. 23)  now  furnish 

T23=T32=T(k)  = 2W'(3  + k2)k  , T33=T(k)k  (0sk<oo)  , (1.26) 

while  all  remaining  components  of  actual  stress  vanish  identically.  We 
shall  henceforth  refer  to  the  graph  of  T(k)  (0sk<oo)  as  the  response 
curve  in  simple  shear.  For  the  special  case  of  the  Neo-Hookean  solid, 
with  the  strain-energy  density  (1.  16),  T(k)=pk  and  this  response  curve 
is  a straight  line.  Bearing  (1.  13)  and  (1.25)  in  mind  one  sees  at  once  that 
the  response  curve  in  simple  shear,  in  turn,  completely  determines  the 
strain-energy  density  for  an  incompressible  elastic  solid  of  the  type 


under  consideration: 
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W(Ij)  = J T(k)dk  (3  £ Ij<oo)  . (1.27) 

0 

Finally,  we  observe  in  connection  with  (1.26)  that  the  Baker- Ericksen 
inequality  (1.  15)  is  equivalent  to  the  obvious  physical  requirement  that 
the  "modulus  of  shear"  2W'(3  + k ) be  positive  at  all  values  of  k. 

When  (1.  15)  is  in  force,  one  confirms  readily  that  the  quasilinear 
second-order  partial  differential  equation  (1.  18)  is  elliptic  at  a solution 
u and  at  a point  (Xj,x^)  if  and  only  if 

W/(3  + k2)  + 2k2W,,(3  + k2)>0  , k=  |vu(xrx2)|  . (1.28) 

On  the  other  hand,  as  is  clear  from  (1.26),  this  inequality  is  satisfied 
if  and  only  if  T7(k)>0,  so  that  the  response  curve  in  simple  shear  has 
a positive  slope  at  an  amount  of  shear  equal  to  the  magnitude  of  the  local 
displacement  gradient  Vufx^x^).  It  follows  that  (1.  18)  cannot  admit  a 
loss  of  ellipticity  if  the  response  curve  in  simple  ohear  represents  a 
strictly  increasing  function*. 

We  conclude  this  section  with  two  examples  of  incompressible 
elastic  solids  within  the  class  of  such  solids  subsumed  in  (1.  13).  The 
first  of  these  concerns  the  subclass  of  "power-law  materials"  introduced 
in  [2],  which  is  governed  by  the  elastic  potential 

W(Il>  = fe{[1  + n(Ir3>]n-  *)  (3£Ij<co)  , (1.29) 

where  |j,  n,  and  b are  positive  material  constants.  In  view  of  (1.26), 


» *If  such  is  the  case,  we  shall  — in  the  present  context  — call  the  elastic 

solid  an  "elliptic  material". 


' 
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the  response  curve  in  simple  shear  here  is  characterized  by 


-.n  - 1 


T(k)  =|_ij"  1 + k (0  £ k<oo)  , 


(1.30) 


being  the  shear  modulus  at  infinitesimal  deformations.  Figure  2 
illustrates  the  response  in  shear  for  various  values  of  the  "hardening 
parameter"  n.  A power-law  material  hardens  or  softens  in  simple  shear 
according  as  n>  1 or  n<l.  The  case  n=l  corresponds  to  the  special 
case  of  the  Neo-Hookean  solid.  For  n^  1/2,  T(k)  is  monotone  strictly 
increasing;  further,  if  n>l/2,  T(k)-»oo  as  k-»oo,  while  T(k)  tends  to  the 
finite  "ultimate  shear  stress"  (iA/2b  as  k-»oo  if  n = 1/2.  On  the  other 
hand,  in  the  range  0<n<l/2,  T(k)  rises  steadily  to  a maximum  at 
k = kQ=V,n/b(l  - 2n)  and  thereafter  decreases  steadily  to  zero  as  k-»oo. 
Thus,  all  solutions  of  the  displacement  equation  of  equilibrium  (1.  18) 
are  bound  to  be  elliptic  everywhere  when  n^  1/2.  In  contrast,  when 
n<  1/2,  a solution  of  (1.  18)  that  involves  displacement  gradients  of  a 
magnitude  exceeding  k^  at  some  point  entails  a local  loss  of  ellipticity. 
The  transition  case  n=  1/2  is  evidently  the  limiting  elliptic  case. 

We  turn  finally  to  a second,  much  more  specialized,  ideal  incom- 
pressible elastic  solid  that  serves  as  a basis  of  the  present  pilot  study. 

The  elastic  potential  defining  this  material  is  given  by 


W(Il)=2(Il-i)  <3sIis4>  ' 
W(I1)=  - ^+2u(Ir3)1/4  (4  £ Ij< oo)  , 


>(1.31) 


in  which  the  positive  constant  m is  again  the  shear  modulus  at  infinitesimal 
deformations.  The  corresponding  response  in  simple  shear  obeys 
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T(k)=uk  (Oikil)  , T(k)=Mk'1/2  ( 1 £ k<oo)  . (1.32) 

Figure  3 displays  the  graph  associated  with  (1.32).  Note  that  the  response 
in  the  present  instance  is  linear  (Neo-Hookean)  up  to  a shear  of  amount 
k=l,  whereas  T(k)  is  monotone  decreasing  for  lsk<oo  and  T(k)-»0 
as  k-»oo.  Clearly,  the  location  of  the  peak  of  the  response  curve  at  k = 1 
is  a matter  of  convenience  only*.  Equation  (1.  18)  at  present  suffers  a 
loss  of  ellipticity  at  a solution  u and  a point  (x^x^)  if  l^ufx^x^)  | > 1. 
Observe  that  while  the  elastic  potential  (1.29)  has  continuous  derivatives 
of  all  orders,  W in  (1.31)  is  merely  once  continuously  differentiable  on 
its  interval  of  definition  [3,  oo);  the  second  derivative  of  this  elastic 
potential  has  a finite  jump  discontinuity  at  Ij  = 4. 

2.  The  small-scale  nonlinear  crack  problem.  Hodograph  transformation. 

Upon  a linearization  with  respect  to  ( Vu | , the  differential  equation 

(1.18)  reduces  to  Laplace's  equation  and  (1.22),  (1.23)  pass  over  into  the 
stress-displacement  relations  appropriate  to  infinitesimal  anti-plane 
shear  deformations.  The  boundary- value  problem  for  u governed  by 

(1.18) ,  (1.20),  (1.21)  under  such  a linearization  becomes  identical  with 
that  characterizing  the  analogous  crack  problem  in  linear  elastostatics. 

This  linear  boundary- value  problem,  in  turn,  is  mathematically  identical 
with  the  one  governing  the  velocity  potential  of  a steady  ir rotational  flow 
of  an  inviscid  incompressible  fluid  past  a flat  plate  of  width  2c  at  a right 
angle  of  attach. 

If  (r,  0)  are  polar  coordinates  at  the  right-hand  crack-tip  (see 

*The  case  in  which  this  peak  occurs  at  k = k^  (k^>0)  is  reducible  to  the 
case  kg=  1 by  means  of  an  obvious  re-scaling. 


f 
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Figure  1)  the  asymptotic  behavior,  as  r-»0,  of  the  displacement  u and  of 

the  shear  stresses  T-  belonging  to  the  well-known  solution  of  the  linearized 

3or 

crack  problem  is  given  by 

* i J/2  • 6 

u~  k(Zcr)  sin  ^ , 

• ° i/->  \- 1 /2  . 8 

T3  i = pu  l ^ikc(2cr)  sin^r  » 

» ° i/o  i-l/2  0 

T32=pu  2~nkc(2cr)  cos  • 

The  normal  stress  vanishes  identically  according  to  the  linear  theory. 

For  our  purposes  it  is  essential  to  remark  that  for  the  Neo-Hookean 
elastic  potential  (1.  16)  the  non-linear  crack  problem  (1.  1 8 ) - ( 1 . 21)  specializes 
exactly  (rather  than  merely  by  linearization)  to  the  linear  problem  of  which 


u is  the  solution.  Further,  for  the  Neo-Hookean  material  coincides 

with  T,  , but  now  no  longer  vanishes.  In  fact,  for  this  special 

3a  33 

material  one  finds  from  (1.23)  and  (2.  1)  that 


T33  = H 


1 2 uk  c 


as  r -*  0 


(2.2) 


According  to  (2.  1),  the  stresses  T3  - and  hence  also  the  correspond- 
ing infinitesimal  shearing  strains  — become  unbounded  as  r-*0.  Consequently 
the  solution  of  the  crack  problem  based  on  linear  elastostatics  cannot  be 
valid  arbitrarily  close  to  the  crack-tips.  It  is  plausible  to  assume,  however, 
that  for  a small  enough  amount  of  shear  k at  infinity,  the  asymptotic 

^See,  for  example,  I.  N.  Sneddon  and  M.  Lowengrub  [ 3 ] , Section  2.6.  Note 
that  the  right  member  in  the  first  of  (2.  1)  is  an  exact  solution  of  Laplace's 
equation  that  has  a vanishing  normal  derivative  along  r>0,  8 = ±it. 


* 
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estimates  (2.  1)  are  realistic  within  a suitably  small  circular  annulus 
centered  at  the  right  crack-tip  (see  Figure  4).  On  the  other  hand,  the 
deformations  and  stresses  in  the  interior  of  the  circular  disk  surrounded 
by  this  annulus  — except  possibly  at  the  tip  of  the  crack,  where  they  may 
not  be  defined  — are  governed  by  the  finite  theory  and  must  match  the 
near  field  predicted  by  (2.  1)  as  the  inner  boundary  of  the  annulus  is 
traversed.  A magnification  of  this  local  state  of  affairs,  which  in  the 
context  of  plasticity  theory*  is  known  as  "small-scale  yielding",  leads 
to  the  following  global  boundary- value  problem  for  the  domain  'S  (see 
Figure  4)  comprising  all  points  exterior  to  the  half-axis  -co<Xj^0, 
x2=0,  which  may  now  be  regarded  as  representing  a crack  of  semi-infinite 
width.  If  (r,  0)  are  at  present  polar  coordinates  at  the  origin,  so  that 

Xj=rcos0  , X£=  r sin  0 (O^rCoo,  , (2.3) 

the  ensuing  "small-scale  nonlinear  crack  problem"  consists  in  finding 
a solution  of  the  partial  differential  equation  (1.  18)  on  &,  subject  to 
the  boundary  conditions 


u 2(x1,0±)  = 0 (— oo  < Xj  < 0) 


(2.4) 


and  the  matching  conditions 


A 


• \ 1 / 2 . 0 
u~k(2cr)  sin  , 

. /->  1/2  .0  . ,,  .-1/2  0 

u i~-kc(2cr)  sin  , u 2~kc(2cr)  cos 


>(2.5) 


J 


1 


See  Rice  [ 4]  . 


' v- : ’ • ..v  sw. . 


as  r -‘oo  (-tt£0<tt)  . 
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Furthermore,  the  desired  solution  u is  to  satisfy  relaxed  regularity 
requirements  strictly  analogous  to  those  imposed  in  conne ction  with  the 
original  nonlinear  crack  problem*. 

When  referred  to  the  dimensionless  variables 


xa=xa/ck2  , r = r/ck2  , 

the  boundar/-value  problem  (1.18),  (2.4), 
non-dimensional  form  in  which  k does  not 


u = u/ck2  , (2.6) 

(2.5)  assumes  the  following 
appear  explicitly: 


[W'(3  + |Vu|2)u  ^=0  on  i 

(2.7) 

u 2(Xj,0±)  = 0 (-ooCXjCO)  ; 

(2.8) 

u~  (2rl  * ^2sin  , 

u j~(2"r)  *^2sin-^  , u ^~(2r)  ^2cos 

9 

2 

M2. 9) 

as  r -*oo  (-it^QCtt)  . 


The  foregoing  small-scale  nonlinear  crack  problem  was  solved 
in  [2]  by  means  of  the  hodograph  method  for  the  limiting  elliptic  case 
n=  1/2  of  a power-law  material  with  the  elastic  potential  given  by  (1.29). 
Reverting  to  the  general  class  of  incompressible  elastic  materials  governed 
by  W = W(Ij),  we  merely  sketch  here  a formal  approach  to  the  problem 
based  on  the  hodograph  method  and  refer  to  [2]  for  a detailed  exposition 


*See  the  regularity  conditions  listed  after  (1.23). 

2 

Subscripts  preceded  by  a comma  are  now  understood  to  indicate  partial 
differentiation  with  respect  to  the  dimensionless  material  coordinates  xQ 
and  V stands  for  the  gradient  operator  with  the  components  d/dxQ  . 
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I I 

of  this  scheme. 

The  underlying  hodograph  transformation  is  the  mapping 

1 

§ =u  (x,,x_)  for  all  (x.,x?)€'#  (2.10) 

from  the  physical  (x^ x2)-plane  into  the  hodograph  plane  spanned  by  the 
rectangular  cartesian  coordinates  (5j/52)*  We  shall  limit  our  attention 
temporarily  to  a material  satisfying  the  ellipticity  condition  T'(k)>0 
(OiSkCoo).  Further,  we  take  for  granted  that  (2.  10)  is  a one-to-one 
continuously  differentiable  mapping  of  'S  onto  its  as  yet  undetermined 
image  A in  the  hodograph  plane.  The  latter  assumption  entitles  us  to 
introduce,  in  place  of  u,  its  Legendre  transform  U as  a new  dependent 
variable: 

U(51,52)  = xpu,p(3c1,x2)-u(x1,x2)  = xp5p-u(x1,x2)  . (2-  11) 


r 


p = ^ = |vu|  . 


(2.  15) 


Under  the  transformations  (2.  10),  (2.  11)  the  nonlinear  boundary- 
value  problem  for  u on  & becomes  a linear  one  for  U on  A,  which  — 
as  shown  in  [2]  — admits  the  explicit  solution 


U(5,,§o)  =pP(p)  coscp  (0< p <oo,  O^cp^ir)  , 


(2.  16) 


where 


p(»)4f  , dt  (o<p<») 

t3W'(3+tZ)  J tZT(t) 


(2.  17) 


On  setting 


Q(p)=-pP'(p)  (0<p<oo)  , 


(2.  18) 


one  infers  from  (2.  16),  (2.12),  in  view  of  (2.  3),  (2.6),  and  (2.  14),  that 


Xj=  r cos  0 = P(p)  - Q(p)  cos  cp  , 


x^=  r sin  0 =-Q(p)  sin  cp  coscp 


(0<p<oo,  O^cp^it)  . 


(2.  19) 


It  can  be  shown  that  for  an  elliptic  material,  equations  (2.  19)  are 
invertible:  for  each  T>0  and  each  0 in  the  interval  [-ir,ir]  there  exist 
unique  values  p=p(r,  0)  and  cp  = cp(r,  0),  with  p>0  and  cp  in  [0,ir], 


The  convergence  of  the  improper  integral  is  assured  by  the  ellipticity 
condition  T '>  0 on  [0,oo). 

2 

A proof  of  this  is  sketched  in  [2]  merely  for  the  limiting  elliptic  power- 
law  material. 


.■'A 


* 


» 

t 


L 
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such  that  (2.  19)  hold.  Once  p(7,  0)  and  cp(r,  0)  are  known,  the  desired 
solution  for  u(xj,x2)  in  the  physical  domain  follows  from  (2.  13)  with 
the  aid  of  (2.  14),  (2.  16),  and  (2.3).  In  this  manner,  bearing  (2.  18)  in 
mind,  one  arrives  at 

u(x1,x2)  = -pQ(p)  coscp  (0<T<oo,  -tt£0<tt)  . (2.20) 

Further,  one  can  verify  a posteriori  that  this  algorithm  indeed  produces 
a solution  of  the  boundary -value  problem  (2.7),  (2.8),  (2.9)  with  the 
requisite  regularity  properties.  In  the  elliptic  case  considered  so  far, 
u so  determined  has  continuous  partial  derivatives  of  all  orders  on  the 
entire  open  region  The  stresses  associated  with  11  are  obtainable 
from  (1.22),  (1.23). 

3 . Elliptic  and  hyperbolic  solutions  of  the  displacement  equation  of 
equilibrium  for  the  special  material. 

Wtih  a view  toward  exploring  the  effect  of  a breakdown  of  ellipticity 
on  the  solution  to  the  small-scale  nonlinear  crack  problem,  we  turn  now  to 
the  special  "pilot  material"  introduced  in  (1.31),  whose  response  to  a 
simple  shear  is  governed  by  (1.32).  When  the  ellipticity  of  the  required 
solution  field  is  not  guaranteed  in  advance,  the  hodograph  mapping  (2,  10) 
may  not  be  one-to-one  and  (2.  19)  need  no  longer  be  uniquely  invertible. 

As  will  be  shown,  however,  even  in  these  circumstances  (2.  16)  to  (2.20) 
may  be  used  to  generate  solutions  of  the  partial  differential  equation  (2.7) 
on  certain  subdomains  of  "S’ . These  solutions,  in  turn,  may  be  combined  to 
construct  a solution  of  restricted  smoothness  to  the  global  boundary- value 
problem  at  hand. 

The  particular  choice  of  the  pilot  material  determined  by  the 
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simple-shear  response  (1.32)  was  motivated  in  part  by  the  fact  that  the 
integral  in  (2.  17)  may  be  evaluated  explicitly  for  this  special  incompressible 
elastic  solid.  Indeed,  (1.32),  (2.  17),  (2.  18)  give 


P(p)  = -^  + | (0<p  si)  , P(p)  = -~  (1*P«»)  , 

2p^  c Jp 


Q(p)  = “V  (0<psl)  , Q(p ) = -p  (lsp<oo)  . 

p vP 


Equations  (2.  19)  may  be  rewritten  as 


x^=  r cos  0 = a(p)  - b(p)  cos  2cp  , 

x^T  sin  0 = -b(p)  sin  2cp  (0<p<oo,  Oscp^ir) 


(3.  1) 


(3.2) 


a(p)  = P(p)--±Q(p)  , b(p)=|Q(p) 


(3.3) 


Thus,  regardless  of  the  specific  choice  of  the  material,  the  points  (x^,*^) 
in  the  physical  plane  that  correspond  to  a constant  positive  value  of  p 
comprise  the  circle 


T(p)  : [xj- a(p)]2+5c|  = b2(p ) , 


(3.4) 


which  is  centered  at  (a(p),0)  and  has  b(p)  as  its  radius. 

Returning  to  the  pilot  material  and  confining  our  attention  to 
0<p<l  for  the  time  being, we  note  that  (3.1)  to  (3.4)  yield  a(p)=3/2, 


b(p)  = l/2p  and 
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'1  2 


2p 


x,  - 4 ~2  cos  2cp  , X£=  - — ^ sin  2cp  (0  < p < 1 , O^cp^ir)  , 


2p 


>(3.5) 


r(p)  : (Xj-  -|)2+x2  = — ^4  (0<P<1) 

4p 


Consequently,  if  (r  j, 0j)  are  polar  coordinates  at  Xj=3/2  , x2=  0 
(see  Figure  5(a)),  so  that 


_ _ 3 _ 

Xj=r  cos  0=  2 + Tj  cos  0j  , 


x^=  r sin  0 = r j sin  0 j (0irj<oo,  -TT^0j<Tr)  , 


one  has  from  (3.5), 


p=(2r1)‘1//2  , cp=-|(0j+Tr) 


Further,  (3.1),  (3.7),  (2.20)  and  0<p<l  now  give 


e, 


_ 1/2  1 

u = e(xj,  x2)  = (2r  j)  sin—  on  £ , 


where  £ is  the  subdomain  of  "2  defined  by 


£ = ( r,  0)  |-^<  TjCoo,  -^<0^^] 


> (3.6) 


(3.7) 


(3.8) 


(3.9) 


and  indicated  in  Figure  5(a).  As  p -»  1 in  the  last  of  (3.5),  T(p)  tends 
to  the  circle 


C : (v|)2+x2=i  °r  ri=  \ ' •fS0i<ff 


(3.  10) 


The  boundary  of  £ thus  consists  of  C and  the  half-line  -oo<x^<l  , 

Since  the  partial  differential  equation  (2.7)  reduces  to  Laplace's 


x2“  °* 


..'■'fit  ; : , 


* 


equation  for  the  pilot  material  in  case  p<  1,  the  displacement  (3.8) 
evidently  satisfies  the  appropriate  specialized  version  of  (2.7)  on  6. 
Moreover,  from  (3 . 8),  (3 . 9)  follows 


|Vu(  = |ve  | =p  = (2r  j)  1 on  6 , 


(3.  11) 


whence  the  solution  (3.8)  is  elliptic  throughout  the  open  region  6. 

Clearly,  |Ve  | -•  1 as  a point  on  C is  approached.  It  is  essential  to 
observe  that  the  elliptic  solution  e is  discontinuous  across  the  half-axis 
-oo<Xj<1  , x^- 0,  rather  than  merely  across  the  semi-infinite  crack 
-oo<Xj<0,  ^2=  O'  1h  fact,  as  is  apparent  from  (2.6),  (2.  1)  and  the 
remarks  following  (2.  1),  the  solution  under  discussion  coincides  with 
the  near  field  of  the  out-of-plane  displacement  for  a Neo-  Hooke  an  material 
and  a crack  tip  situated  at  7 j = 0,  rather  than  at  the  origin  r = 0. 

We  proceed  now  to  the  case  p>  1 . In  this  instance  (3 . 1)  to  (3.4) 
furnish  a(p)  = 3/2Vp  , b(p)=  1/2,/p 


x.  = r cos  0 = — —(3  - cos  2cp)  , 

1 Vp 


x.7-  r sin  0 = - — — sin  2cp  (l<p<oo,  0 s cp ^ it)  , 


(3.  12) 


r(p)  :(xr7^)  +xz  = ~k  (1<p<oo)  • 


Upon  referring  the  last  of  (3.  12)  to  the  polar  coordinates  (7,  0)  and 
solving  the  resulting  quadratic  equation  in  Vp  » °ne  is  led  to 


Jp  = ”[3  cos  0±R(0)]  , 

27 


(3. 13) 


provided 


R(9)  = - 9 cos  0-8 


(3.  14) 


On  the  other  hand,  (3.  13),  (3.  14)  and  the  first  two  of  (3.  12),  after  some 
algebraic  manipulations,  are  found  to  imply 


-2 Jl  sin  6 

co  s cp  — o 1/2 

[ 1 + 3 sin  0±R(9)cos0] 


cos  0 ± R(9) 

sm  cp  = 7 4 f72 

72[l + 3 sin  9±R(9)  cos  0]  ' 


(3. 15) 


Clearly,  both  values  of  «/p  supplied  by  (3.  13)  are  real  and  positive 


if  and  only  if 


cos  9 1*7879  or  -0o*6*9q,  0q=cos  l(j8/3)  . (3.16) 


Further,  since  p>l  at  present,  (3.  13)  demands  that 


3 cos  0±R(0)>  2r  . 

Now  let  M and  Q be  the  two  subdomains  of  '3'  defined  by 

V = {(r,  9)|0<r<|[3  cos  0 + R(0)]  , Os  | 0 | < 0Q}  , 

V = [(r,  0)|O<r<-|[3cos0-R(0)]  , 0<  | 0 | < 0O}  , 


(3.  17) 


(3.  18) 


which  are  shown  in  Figure  5(b)  and  Figure  5(c),  respectively,  and  which 
evidently  conform  to  (3.  16),  as  well  as  to  the  inequality  (3.  17)  with  the 


^The  reason  for  excluding  points  on  the  ray  9 = 0 from  V will  emerge 
pre  sently . 


£5  «. 
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appropriate  sign  alternative. 

For  the  purpose  of  examining  the  geometry  of  the  two  open  point  sets 

V and  V we  note  first  that  T(p)  in  the  last  of  (3.  12),  as  p-*l,  tends 

+ 

to  circle  C given  by  (3.  10).  Next,  let  C and  C be  the  two  complementary 
circular  arcs  of  C defined  by 


c = {(Fj, e1)|r1=  j , OS  | Q j | ^BQ+fl  * 

c = {(?1,e1)|71=| , e0+|<  1 9j | =stt}  , 


M3.  19) 


j 


and  observe  on  the  basis  of  (3.6),  (3.  14),  (3.  16)  that  these  two  arcs  may 
also  be  described  by 


t = {(r,  0)|r=  -|[3  cos  0 + R(0)]  , 0 ^ | 0 | ^ 0Q}  , 

C = f(r,0)|r  = -|[3cos0-R(0)]  , 0 ^ | 0 | < ©q}  . 


> (3 . 20) 


+ 

It  is  clear  from  (3.  18),  (3.20)  that  the  boundary  of  V consists  of  the  two 
ray-segments  0sr^^(  9 = ±00  together  with  the  circular  arc  C,  as 

indicated  in  Figure  5(b);  in  contrast,  the  boundary  of  V — displayed  in 
Figure  5(c)  — is  now  seen  to  be  composed  of  the  same  two  ray- segments, 
the  circular  arc  C,  as  well  as  the  segment  O^T^l,  0 = 0 of  the 


x j - axis . 

Recalling  from  (3.  1)  that  Q(p)=  l/2Vp  when  p>l  and  substituting 
from  (3.  13)  and  the  first  of  (3.  15)  into  (2.20),  we  arrive  at  the  two  out- 
of-plane  displacement  fields: 
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u = h(x.  ,x?)  = UcosJ  + RMUl^inJ 

v/2  rtv-oi  e + R(0)]  ' 


+ 

on 


» = h(x,  ,s,)  ° Jas--9~  R<9>1 3 /2? 

-s/2  r[  cos  0 - R ( 0) ] 


on  V 


(3.21) 


with  R(0)  given  by  (3.  14).  Both  h and  h are  real-valued  functions  that 


have  continuous  partial  derivatives  of  all  orders  on  their  respective  domains 


of  definition;  further,  each  satisfies  the  displacement-equation  of  equili- 


brium (2.7),  specialized  for  the  pilot  material  and  p>l,  as  is  readily 


verified  by  substitution.  Also,  a direct  computation  confirms  — in  agree- 


ment with  (2.  15)  and  (3.  13)  — that 


|Vu|  = |vh|  = -L[3cos  0 + R(6)]2>  1 on  l , 


|vu|  = |vh|  = — i-[3cos  0 - R(0)]2>  1 on  V , 


(3.22) 


so  that  h and  h are  both  hyperbolic  solutions  throughout  the  open  sets 

+ - ._+ . 

V and  W , respectively.  From  (3.22)  and  (3.20)  one  infers  that  |vh|  -*  1 

+ + 

as  a point  on  the  boundary- segment  C of  V is  approached;  similarly. 


|Vh|  -*  1 at  every  point  of  the  circular  arc  C. 


Finally,  (3.  14)  and  the  second  of  (3. 21)  reveal  that 


h(xj,  0±)  = ± =-  (0<xj<l) 


(3.23) 


Consequently,  the  second  of  the  foregoing  two  hyperbolic  solutions  has  a 


finite  jump  discontinuity  across  the 


jnt  0<Xj<1  of  the  x^-axis. 


..'F!  •*-'&<*  * 


, 
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4.  Construction  of  a solution  to  the  small-scale  nonlinear  crack  problem 
for  the  pilot  material. 

In  this  section  we  attempt  to  piece  together  the  solutions  of  the 
partial  differential  equation  (2.7)  established  in  Section  3 in  order  to  con- 
struct a solution  — continuous  and  piecewise  smooth  on  ^ — of  the  small- 
scale  nonlinear  crack  problem  (2.7),  (2.8),  (2.9)  for  the  pilot  material. 

To  this  end  we  note  first  that  we  have  available  only  the  elliptic 
solution  e on  the  complement  of  W with  respect  to  the  cut  plane  U and 
only  the  hyperbolic  solution  h in  the  interior  of  the  circle  C (see  Figure  5). 

Further,  as  is  easily  seen  from  (3.8),  (3.21)  and  (3.6),  (3.14),  both  e 

— + _+  + 
and  Ve  join  continuously  with  h and  Vh  across  the  circular  arc  C 

+ 

separating  the  elliptic  domain  t and  the  hyperbolic  domain  V.  Thus  e 
and  h supply  a solution  of  the  requisite  smoothness  to  the  differential 
equation  (2.7)  at  all  points  of  f?  that  lie  outside  the  boundary  of  K . 

Since  h is  the  only  one  of  the  three  solutions  available  on  27  that 
is  continuous  across  the  segment  0<Xj<  1 of  the  x^-axis,  it  follows 
that  the  solution  to  the  boundary- value  problem  (2.7),  (2.8),  (2.9)  we 
are  trying  to  assemble  must  coincide  with  h in  the  vicinity  of  this  seg- 
ment. On  the  other  hand,  the  transition  from  h to  e across  the  ray- 
segments  0^Ts^2,  0 = ±0O  is  readily  found  to  involve  a jump  discontinuity 

in  u. 

The  foregoing  considerations  suggest  that  we  dismiss  h altogether 

+ 

and  seek  to  "match"  h and  e across  two  as  yet  undetermined  curves 
J and  V'  (Figure  6)  that  are  symmetrically  situated  with  respect  to  the 
Xj-axis,  start  at  the  origin,  terminate  at  the  points  r = , 9=±9q»  and 

otherwise  lie  wholly  in  27.  This  matching,  moreover,  must  assure  the 


continuity  of  u and  of  the  traction  across  the  two  curves  to  be  found.  As 
will  emerge  shortly,  there  exist  unique  "shocks"  (shock-lines)  J and 
with  the  desired  properties;  as  either  of  these  lines  is  traversed  u and  its 
tangential  derivative  remain  continuous,  while  the  normal  derivative  of 
u suffers  a finite  jump  discontinuity  — despite  the  prevailing  traction  con- 
tinuity across  J and 

The  final  solution  of  the  small-scale  yielding  problem  thus  obtained 

+ 

coincides  with  the  hyperbolic  solution  h on  the  finite  domain  ft'.,.  whose 
boundary  is  the  closed  curve  composed  of  J , J',  and  C (see  Figure  6); 
it  coincides  with  the  elliptic  solution  e on  the  domain  t ^ consisting  of 
all  points  of  £ exterior  to  the  boundary  of  ft^  . Since  'Fj~‘r  and  0j~9 
as  T-*co,  it  is  clear  from  (3.8)  that  the  composite  solution  just  described 
satisfies  the  conditions  (2.9)  at  infinity,  in  addition  to  the  boundary  con- 
ditions (2.8),  which  are  obviously  met.  While  the  solution  here  alluded  to 
clearly  possesses  the  required  (relaxed)  smoothness  on  &,  we  shall  have 
to  show  eventually  that  u so  determined  remains  bounded  as  "r-»0. 

We  proceed  now  to  the  determination  of  the  upper  shock  J , from 
which  of'  follows  by  symmetry.  Since  the  hyperbolic  solution  h has  been 
discarded,  we  shall  henceforth  write  h in  place  of  h.  According  to  (3.8) 
and  (3.21),  the  two  solutions  of  (2.7)  to  be  matched  across  the  desired 


curve  a/  are: 
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where 


f(e)  = J-t3cos9  + R(9)]3/2sin9  # R{Q)  J9cosZ  q.  8 . (4.2) 

72  [cos9  + R(9)J 


The  continuity  of  the  displacement  u across  J therefore  demands  that 


1 2 r 


on 


(4.3) 


in  which  r^  and  9j,  on  account  of  (3.6),  may  be  regarded  as  functions 
of  "r  and  9.  We  shall  show  that  (4.3)  actually  admits  a unique  solution 

^:7=r(9)  (0  £ 9 £ 9Q)  , r (0)  = 0,  r(9Q)  = Jz  (4.4) 

and  that  the  emerging  shock  J meets  all  the  conditions  stated  earlier. 

With  a view  toward  proving  this  claim,  we  temporarily  take  for 
granted  the  existence  of  a solution  (4.4)  to  (4.3)  and  infer,  with  the  aid 
of  (3.6),  that  then  necessarily 

r^sin29  - 2r3f2(9)cos  9 + 3r 2f2(9)  - f^(9)  = 0 . (4.5) 

Next,  we  introduce  the  auxiliary  function 

C = C(0)  = Scos^ +^R(9)  (4-6) 

and  use  (4.2)  to  conclude,  after  lengthy  elementary  manipulations,  that 
(4.  5)  may  now  be  written  as 

(C-  |)2{c4+C3+|c2-(4C+1,7 — -2}  = 0 . (4.7) 

L L 4 [cos9  + R(9)]  j 

But  (4.  6)  and  (3. 20)  imply  that  the  first  factor  in  (4.  7)  vanishes  if  and  only 
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if  r = r(9)  and  0 are  the  polar  coordinates  of  a point  on  the  circular 

arc  £ and  hence  this  factor  cannot  vanish  on  J.  The  presence  of  the 

+ 

root  Q = 1/2  of  (4.7)  merely  reflects  the  continuity  of  u across  C 
mentioned  previously.  It  follows  that  the  function  within  braces  in  (4.7) 
must  vanish  on  so  that 


0(0  = 


c4+c3+|c2 


. zn 

sin  9 


^ + * [ cos  0 + R( 9) ] ^ 


M(e) 


(0<;  0 5=  0O)  on  J . 


>(4.8) 


■ 

The  function  ill  appearing  in  (4.8)  is  strictly  increasing  on  [0,  0Q],  with 
\|r(0)=0  and  \|i(0q)  = 1/8,  whereas  the  rational  function  0 on  the  left 
side  of  this  equation  is  easily  seen  to  increase  steadily  from  0(0)  = 0 
to  0(1/2)  = 1/8  as  Q ranges  over  the  interval  [0,  1/2].  Consequently, 
for  each  given  0 in  [O,0Q]  there  exists  precisely  one  root  £ = £(0)  of 
(4.8).  If  we  now  employ  (4.6)  to  define  r(0)  and  reverse  the  steps  that 
led  from  (4.3)  to  (4.8),  we  infer  the  existence  and  uniqueness  of  the  shock 
J described  by  (4.4). 

Although  (4.  8)  has  not  been  solved  explicitly  for  Q = Q(Q),  the 
numerical  calculation  of  r(0)  (Ofi0S  0Q)  from  (4.8),  (4.6)  presents  no 


difficulties.  The  results  of  this  computation  are  plotted  in  Figure  6. 
Further,  equations  (4.8),  (4.6)  enable  one  to  extract  the  asymptotic 
behavior  of  the  shock  J near  its  endpoints.  One  finds  in  this  manner  that 
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As  is  easily  verified,  (4.9)  assures  that  J is  tangential  to  the  Xj-axis 
at  the  origin  and  tangential  to  the  circle  C at  r = */2  , 0 = 0q- 

Our  next  objective  is  to  establish  the  continuity  of  the  tractions 
across  the  shock  J . In  this  connection  we  rely  on  (1.24),  (1.  19),  and 
(2.  6)  to  see  that  the  only  non- vanishing  nominal  traction  component  along 
J is  given  by 

s = s3=2W/(3+  |Vu|2)||  , (4.10) 


where  u = e on  the  elliptic  side  of  J,  while  u=h  on  the  hyperbolic  side. 
Denoting  the  scalar  traction  s on  the  two  sides  of  J by  sg  and  s^, 
respectively,  we  draw  from  (4.  10)  and  (1.26),  (1.32)  that 


3e  l—,  | on 

se  = ^9K  ' sv,  = MVhl  ^7  • 


■3/2 


3h 

9n 


(4. 11) 


Accordingly,  the  tractions  are  continuous  across  J if  and  only  if 


sh-se=0  or  { |vhr3/2Vh-Ve|.n  = 0 on  J , (4.12) 

provided  n is  a unit  normal  vector  of  J and  the  dot  signifies  scalar 
multiplication.  Keeping  in  mind  that  the  displacement- continuity  condition 
(4.3),  which  determined  the  shock  J,  asserts 


one  has 


h - e = 0 on  J , 


Vh-  Ve 
|vh-  Ve 


on  J . 


(4.  13) 


(4.  14) 


The  traction- continuity  condition  (4.  12)  is  therefore  equivalent  to 


(4.  15 


|vh|  1/2+  |ve  |2- {l  + |vh|"3/2|vh*ve  = 0 on  J . 

Aiming  at  a verification  of  (4.  15),  we  first  record  the  polar  com- 
ponents of  Vh,  which  follow  from  the  second  of  (4.  1): 


_ — f(Q)  — i-ll  - — /(fl) 

9r  -21'9'  ’ 7 90  " — 2*  9 


(4.  16) 


Similarly,  the  first  of  (4.  1),  in  conjunction  with  (3.6),  enables  one  to 
arrive  at 


9e 

90 


J2*l 


cos 


Equations  (4.  16),  (4.  17)  reduce  (4.  15)  to 


(4.  17) 


=[f2(0)  + f2(0)]l/4 


with 


- -{l  + r3[f2(9)  + f2(6)r3/4j-F(e)  = 0 on  J , (4.18) 

r ^/2¥\^ 


F(9)  = /(0)cos(e--y-)-f(e)sin(e--y-)  . (4.19) 


Now,  (3.6)  and  (4.3)  justify 


•=[r: 


f (0)sin  0j  -,1/3 


2sin  0sin  (0^/2) 


:r  [ 


f(0)sin  0 


«/2sin(8j  /2)sin  0^ 


n2/3 

J 


on  J . (4.  20) 


Here  and  in  what  follows  l stands  for  the  derivative  of  f. 


I 
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Finally,  substitution  from  (4.  20)  into  (4.  18)  and  subsequent  use  of  the 
definition  (4.2)  of  f ( 0) , after  a lengthy  computation,  confirm  that  (4.  18) 
is  indeed  satisfied.  This  completes  the  verification  of  the  traction-con- 
tinuity requirement  (4.  12). 

We  have  yet  to  show  that  the  solution  u constructed  in  this  section 
remains  bounded  as  T-»0.  According  to  the  first  of  (4.  1),  the  restriction 
of  e to  the  final  elliptic  domain  (Figure  6)  clearly  meets  this  condition. 
It  therefore  suffices  to  prove  that  the  hyperbolic  solution  h is  bounded  on 
the  point-set  9 consisting  of  all  points  (r,  Q)  with  Y>0,  0^0  that  belong 
either  to  V or  to  J . Note  that  9 is  the  closure  of  the  upper  half  of 
the  final  hyperbolic  domain  , with  the  exception  of  the  origin. 

Since  J is  tangential  to  the  x^-axis  at  the  origin*,  it  follows  that 
every  ray  from  this  point  to  a point  in  9 intersects  J precisely  once, 
whence 

T^r(0)  for  all  fr,  9)  in  9 . 

We  recall  next  that  £(9)  in  (4.8)  has  the  range  [0,  1/2],  so  that 

co/e^ei  <oseiV 

and  (4.6)  yields  the  lower  bound 

g(9)  a /|  [ 3 cos  9 + R.(9)  J sin  e (Os0^0  ) (4.23) 

73  [ cos  9 + R(9 )]  0 

But  (4.  21)  and  the  second  of  (4.  1)  give 
*Refer  to  the  first  of  (4.9). 


(4.21) 


(4.22) 
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i 


h(*,.x2)  = ^^  on  9 , (4.24) 

while  this  inequality,  together  with  (4.23)  and  (4.  2), leads  to 

h(xj,  ^ fa  3 cos  9 + R(9)][cos  9 + R(0)]}1/2^  Jl  on  9 . (4.25) 

This  establishes  the  required  boundedness  of  h,  and  thus  of  u , as  r-*0. 

In  contrast,  tfu  fails  to  remain  bounded  as  the  origin  is  approached 
from  within  K, , as  is  apparent  from  the  first  of  (3.22),  which  gives 

| Vu(Xj,  0)  | ~ as  Xj-*0+  . (4.26) 

X1 

One  shows  easily  that  |vu|  -oo  as  r-0  along  any;  regular  arc  lying  in  V* 
or  along  J.  In  addition,  it  should  be  pointed  out  that  (4.  1)  furnishes 

TT(x j , 0i)  = ±*JZr ^ (-oo<x1<0)  , u(xlf0)  = 0 (OCXjCco)  , (4.27) 

which  reveal  a displacement  discontinuity  at  the  origin. 

There  is  one  further  requirement  that  has  been  proposed  for 
elastostatic  fields  with  discontinuous  deformation  gradients  which  has  not 
been  mentioned  so  far  in  the  present  paper.  This  condition  insures  the 
dissipative  character  of  the  shocks  inherent  in  such  a field  and  is  an 
analogue  of  the  entropy  inequality  for  gas-dynamical  shocks. 

The  elastostatic  inequality  to  which  we  are  alluding  was  first 
deduced  in  [5]  with  limitation  to  rectilinear  shocks  in  piecewise  homogeneous 
plane  deformations;  it  was  later  extended  to  curved  shocks  and  general 
non-homogeneous  deformations  in  [6],  where  the  special  form  assumed  by 
* this  condition  in  the  particular  case  of  anti-plane  shear  was  given  separate 

» 

I I 

...  . ).  i 
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attention*.  A thermodynamic  argument  in  support  of  the  dissipation 
condition  for  compressible  elastic  materials  is  produced  at  the  end  of  [6]. 

The  application  to  the  solution  field  under  consideration  of  the 
dissipation  inequality  referred  to  above  is  contingent  upon  a quasi- static 
interpretation  of  this  field,  in  which  the  amount  of  shear  k applied  at 
infinity  plays  the  role  of  "time"  — or,  more  appropriately,  of  a history 
parameter  . In  this  connection  the  corresponding  velocity  V(x,  k)  of 

a point  x on  the  instantaneous  shock  «^(k)  relative  to  the  undeformed 

3 

configuration  is  of  importance  . 

In  view  of  (4.4)  and  (2.6),  the  physical  position  vector  of  a point 
on  »^(k)  is  given  by 

x = x(e,k)  = ck2r(0)a(0)  (O^0^0J  , (4.28) 

/v  rw  (J 

if  a(0)  is  the  radial  unit  base  vector  associated  with  the  polar  coordinates 
(r,  0 ).  Hence 

V(x,k)=|^x(0,k)  = |x(0,k)  (OS0^0O)  . (4.29) 

It  follows  from  (4.28),  (4.29)  and  properties  of  the  function  ? established 
earlier  that  V(x,  k)  always  points  into  the  hyperbolic  domain  , so  that — 
as  k increases  — material  points  of  the  body  traverse  the  shock  »/(k)  from 
the  hyperbolic  "upstream"  to  the  elliptic  "downstream"  side  of  »^(k).  The 


*See  Section  6 of  [6]. 

2 

Thus  as  "time"  k increases,  the  body  undergoes  loading,  rather  than 
unloading. 


3 — 

By  virtue  of  the  symmetry  of  the  solution  u 

following  considerations  to  the  upper  shock. 


at  hand  we  may  confine  the 


r 


relevant  dissipation  inequality  now  takes  the  form  of  the  jump  condition 


[H]+*0  on  J , 


(4.30) 


in  which  the  function  within  brackets  is  defined  by 


H = - W(3  +|vu|2)  + 2W/(3+|vu|  2)(|£) 


(4.  3 1) 


and  [h]*  stands  for  the  difference  between  the  limiting  value  of  H on 
the  upstream  (or  positive)  side  of  J and  its  limiting  value  on  the  downstream 
(or  negative)  side.  On  account  of  (4.31),  (4.  10),  (2.6),  and  the  continuity 
of  the  tractions  across  J,  we  arrive  at 

[H]t=-[w]t+s[ff]  on  J , W=  W(3+  |v5| 2)  . (4.32) 

Recalling  that  u = h on  the  upstream  side,  while  u = e on  the 
downstream  side  of  J,  we  now  express  the  right-hand  member  in  the  first 
of  (4.32)  in  terms  of  vh  and  Ve.  To  this  end  we  first  invoke  (1.31), 

(1.  19),  and  (2.6),  which  give 


[W]l=^{4|vh|  1/2-3-lve(2]  , 


(4.33) 


whereas  (4.  11)  and  the  first  of  (4.  12)  enable  us  to  write 


3f = 0 I Vb | " 3 7 2 vh-n(vh*n-  Ve-n) 


(4.34) 


Furthermore,  another  appeal  to  (4.  11),  (4.  12)  yields 


^See  (4.  19)and  (6.  13)  in  [6] . 

2 

Here  9u/0n  = Vu*n  and  n is  a unit  normal  vector  of  J. 


V *-.V;  ' : <*-  . \ 
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Ve*  n = |vhr3/2Vh*n  . 


(4.35) 


Using  this  result  in  (4.34)  and  thereafter  substituting  from  (4.33),  (4.34) 
into  (4.32),  we  are  finally  led  to 

[H]t=^{lve|2+3-4lvh|1/2+  ^\^2[l-  |vh|-3/2]j  . (4.36) 


We  turn  next  to  a verification  of  the  dissipation  inequality  (4.30) 
at  the  endpoints  of  the  shock  J.  Since  | Vh | > 1 , (4.36)  implies  that 


[H]^  |ve  |2-  2 | Vh | 1/2+3}  on  J . 


(4.37) 


But  according  to  (3  . 1 1)  and  (3 . 22),  |Ve  | is  bounded  and  |vh|-oo  as 
T-*0.  Consequently,  (4.30)  is  certain  to  hold  along  J in  a neighborhood 
of  the  origin. 

The  confirmation  of  (4.30)  near  the  right  endpoint  of  J , at  "r=,/2  , 
9 = 9q  , unfortunately  turns  out  to  be  a very  cumbersome  task,  which 
requires  a considerable  improvement  of  the  second  estimate  in  (4.9).  The 
requisite  refinement  of  this  asymptotic  result  is  found  to  be 


f(0)  =72  - 23/*(0o-  0)1/2+  |(0Q-  9)  - 226x2"  "'*(0O-  0)" 


+ O{(0Q-  0 ) } as  0-0o  , 


(4.38) 


and  calculations  far  too  elaborate  to  be  included  here  eventually  yield 


[H]!=-33x21/4(0o-0)3/2+o{(0o-0)2}  as  0-0Q  . 


(4.39) 


’ ■ 
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It  is  now  evident  that  (4.30)  also  holds  in  a neighborhood  of  the  right 
endpoint  of  J.  Numerical  calculations,  which  were  carried  out  on  the 
basis  of  (4.36),  reveal  that  the  dissipation  inequality  is  in  fact  satisfied 
at  all  points  of  the  shock. 

In  defining  the  pilot  material  through  (1.31)  we  assumed  for  con- 
venience that  the  peak  of  the  response  curve  in  simple  shear  (Figure  3) 
is  situated  at  k=  1.  If  this  normalization  is  abandoned  and  the  maximum 
of  T(k)  is  permitted  to  occur  instead  at  k = k^  (k^X)),  (1.32)  give  way  to 

T(k)=nk  (0£k£k0),  T(k)  =^ik0(k/k0)"  (kQ^k<oo)  , (4.40) 

while  (1.31)  undergo  an  analogous  modification  in  accordance  with  (1.27). 

It  is  easily  seen  that  in  these  circumstances  the  physical  problem  governed 
by  (1.18),  (2.4),  (2.5)  again  reduces  to  the  dimensionless  version  (2.7), 
(2.8),  (2.9)  provided  the  scaling  (2.6)  is  replaced  by 

xa=  k^x^/ck2  , r = k^r/ck2  , u = kQu/ck2  . (4.41) 

Since  the  solution  to  the  small-scale  nonlinear  crack  problem 
deduced  here  involves  shocks,  along  which  the  displacement  gradients 
suffer  jump  discontinuities,  it  is  natural  to  wonder  whether  it  is  possible 
to  prove  the  non-existence  of  a smoother  solution  to  this  problem.  We 
have  not  been  able  to  supply  such  a proof;  nor  have  we  dealth  with  questions 
related  to  the  stability  of  the  solution  found. 

5.  Numerical  results  for  the  displacements  and  stresses.  Generalizations 
of  the  constitutive  assumption.  Discussion. 
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the  displacement  field  of  the  solution  to  the  pilot  problem  deduced  there 
is  supplied  by 


0 

u -JZr  j sin  y on  6^  , u = on  , 


(5.  1) 


where  f (9 ) is  given  by  (4.2).  Here  6.,.  and  V,,.  are  the  final  elliptic 
and  hyperbolic  domains  defined  in  Section  4 and  depicted  in  Figure  6. 
Equations  (5.  1)  together  with  the  relations  (3.6)  between  the  two  pairs  of 
polar  coordinates  (r^,9j)  and  (r,  9)  maybe  used  to  calculate  u(r,9). 
Figure  7 shows  curves  displaying  the  angular  dependence  of  u(r,  9)  for 
four  fixed  values  of  "r.  Note  that  the  kink  in  the  curves  corresponding  to 
"r  = 1/4  and  "r  =^2/2  is  due  to  the  shock.  Since  the  circles  "r  =7/4  and 
r = 9/4  do  not  intersect  the  shock,  the  remaining  two  curves  in  Figure  7 
are  free  of  any  discontinuity  in  slope. 

From  (5.  1)  follows 


0<  9 ^ TT 


lim  u(r,  9) 
T— 0 


0 for  9 = 0 
-J3  for  -tt^9<0  , 


(5.2) 


which  accounts  for  the  dashed  line  in  Figure  7.  Thus,  in  first  approximation, 
the  material  above  and  below  the  crack,  near  its  tip,  translates  rigidly 
in  opposite  directions  (parallel  to  the  plane  of  the  crack).  The  displacement 
discontinuity  at  the  crack  tip  inherent  in  (5.2)  was  mentioned  earlier  in 
connection  with  (4.  27). 


We  proceed  now  to  the  determination  and  discussion  of  the  stresses 
associated  with  the  solution  under  consideration.  In  view  of  (1.  19),  (1.  22), 
(1.23),  the  non-trivial  components  of  the  prevailing  actual  stress  field, 
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referred  to  the  cylindrical  coordinates  (r,  8,  z)  = (r,  6,  x^),  are: 


Tr«=2Wl'Il>f?'  T6.=  2Wl>IX>?ff  ' %z=2W'«Il>l7“l 


//t  vA 


ij=3+ ivur  , 


> (5.3) 


where  W is  the  elastic  potential  of  the  pilot  material  defined  in  (1.31). 
A direct  computation  relying  on  the  scaling  (2.6),  the  dimensionless 
displacements  (5.1),  and  involving  the  polar -coordinate  relations  (3.6), 
ultimately  yields  the  results: 


T 

zz 


-_l i- 
271 


on  e*  , 


>(5.4) 


which  — together  with  (3.6)  — determine  the  cylindrical  components  of 
the  Cauchy  stress  tensor  on  the  elliptic  domain  as  functions  of  r 
and  0 alone;  similarly. 


Tr2  = -Urf(e)[f2(0)  + f2(0)]"3/4  , 

T0z=pr/(0)[f2(0)  + /2(0)]"3^4  , 

t =£[f2(0)  + f2(0)]1/4  on  V*  , 
zz  r * 


>(5.5) 


J 


with  f(0)  furnished  by  (4.2). 


Figures  8,9,  10  contain  graphs  depicting  the  dependence  of  the  stresses 
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t , Tfl  , t upon  the  polar  angle  0 for  the  same  four  fixed  values 
of  7 employed  in  Figure  7,  except  that  the  results  for  7=9/4  were 
omitted  in  Figure  10  in  order  to  avoid  an  excessively  crowded  diagram. 
The  curves  appropriate  to  7 = 1/4  and  7 =J2/Z  in  the  figures  under 
discussion  display  finite  jump  discontinuities  at  values  of  0 that  are 
angular  polar  coordinates  of  the  point  of  intersection  of  the  corresponding 
circle  7 = constant  with  the  shock.  In  contrast,  the  graphs  based  on 
7 = 7/4  reveal  merely  a discontinuity  in  slope  at  the  value  of  0 that 
corresponds  to  the  intersection  between  the  circles  C and  7=7/4. 

At  this  point  of  intersection  |vu|  =1  and  the  singular  behavior  just  men- 
tioned can  be  traced  to  the  kink  in  the  shear- response  curve  (Figure  3) 
at  k=  1.  This  kink  is  also  responsible  for  a discontinuity  in  the  curvature 
of  the  displa''  ment  curve  for  7=  7/4  in  Figure  7 at  the  particular  value 
of  0 referred  to  above.  Finally,  the  curves  associated  with  7 = 9/4 
in  Figures  8,9  are  entirely  smooth  since  this  circle  lies  wholly  in  the 
domain  of  ellipticity  of  the  solution. 

If  the  maximum  of  T(k)  occurs  at  an  amount  of  shear  k^  , so  that 

(4.40)  hold  and  the  scaling  (4.41)  is  in  force,  the  graphs  in  Figures  8,9 

represent  t /pk  , and  tq  /pk  , while  those  in  Figure  10  become 
rz  u tjz  u 

2 

plots  of  Tzz/pkQ  . 

Equations  (5.4),  (5.5)  imply  that  the  shear  stresses  Trz>  tqz 
remain  bounded  as  7-»0  and  tend  to  zero  as  the  crack-tip  is  approached 


from  within  the  hyperbolic  domain  . This  conclusion  is  consistent 
with  (4.26)  and  the  fact  that  T(k)-»0  as  k-«oo  in  a simple  shear  defor- 
mation of  the  special  material  under  consideration.  On  the  other  hand. 
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the  axial  normal  stress  t , though  bounded  in  any  elliptic  neighborhood 

z z 

of  the  origin  7=0,  becomes  unbounded  as  7-»0  from  within  . This 

singular  behavior  of  r was  to  be  anticipated  in  view  of  the  last  of  (1. 26), 

according  to  which  the  normal  stress  induced  by  a simple  shear 

remains  finite  as  the  amount  of  shear  k grows  beyond  bounds  only  if 

1/2 

kT(k)  tends  to  a finite  value  as  k-»oo,  whereas  kT(k)~|ik  in  this 

limit  for  the  pilot  material*.  The  preceding  remark  suggests  that  the 

singularity  of  t at  the  crack-tip  would  not  arise  if  the  second  of  the  re- 
z z 

lations  (1.32)  defining  the  material  were  suitably  modified.  Unfortunately, 
this  objective  is  incompatible  with  the  hodograph  scheme  we  have  used  to 
deduce  the  solution  to  the  problem  at  hand  since  the  improper  integral  in 
(2.  17)  is  no  longer  convergent  when  kT(k)  = 0(l)  as  k-»oo. 

Keeping  in  mind  that  |Vu|<l  on  and  | V u | > 1 on  one 

easily  infers  from  (1.19),  (1.23),  and  (1.31)  that 

W = 4t  on  e„.  , W=2t  - on  K.  . (5.6) 

The  behavior  of  the  strain-energy  density  associated  with  the  solution 
being  discussed  is  therefore  immediate  from  that  of  the  axial  stress  t 

z z 

In  particular,  W remains  bounded  at  the  crack-tip  except  when  the  latter 
is  approached  from  within  . 

An  analysis  strictly  parallel  to  that  summarized  in  the  foregoing 
sections  leads  to  a solution  of  the  analogous  small-scale  nonlinear  crack 
problem  for  an  incompressible  elastic  material  (of  the  type  characterized 
by  (1.  13))  whose  response  in  simple  shear  obeys 

*See  the  second  of  (1.32). 
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T(k)=|jk  (Osk^l),  T(k)=(ik"°  (Uk<oo)  , (5.7) 

in  which  a is  a material  constant  restricted  by 

0 £ a<  1 . (5 . 8) 1 

The  constitutive  relation  (5.7)  is  evidently  a generalization  of  (1.32)  and 
reduces  to  (1.32)  for  a=l/2.  The  special  case  a = 0 corresponds  to  a 
shear  response  in  which  the  shear  stress  T(k)  rises  linearly  to  a value 
|j  at  k = 1,  and  remains  constant  at  this  value  as  k increases  from  unity. 

For  materials  governed  by  (5.7),  (5.  8), the  differential  equation 
(1.  18)  again  suffers  a loss  of  ellipticity  at  a solution  u and  at  a point 
(xj,x2),  provided  |vu(Xj,  x2)  | > 1. 

On  applying  the  hodograph  scheme  to  the  small-scale  nonlinear 
crack  problem  for  a material  with  the  shear  response  (5.7),  (5.8),  one  is 
once  again  led  to  elliptic  and  hyperbolic  solutions  of  (1.  18)  on  certain 
subdomains  of  7}  . In  fact.  Figure  5 remains  applicable  in  the  present 
circumstances,  provided  the  circle  C is  replaced  by  one  of  radius  1/2 
centered  at  Xj=  |(l  + a)/(l  - a)  , x2=  0 , while  the  angle  0Q  is  now 
determined  by  sin  6Q=  (1  - »)/(!  + <*).  Further,  the  formulas  (3 . 8)  and 
(3.21)  for  the  elliptic  and  the  two  hyperbolic  solutions  at  present  undergo 
modifications  which  will  not  be  spelled  out  explicitly.  As  in  the  special 
case  a = 1/2,  the  elliptic  solution  can  be  matched  to  the  appropriate 
hyperbolic  solution  across  two  symmetrically  situated  shocks,  if 
0<a<l.  Thus  Figure  6,  which  illustrates  the  elliptic  and  hyperbolic 
domains  of  the  final  solution  for  or  = 1/2,  has  a qualitatively  similar 

*The  restriction  cr<  1 assures  the  convergence  of  the  integrals  in  (2.  17). 
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[ 


counterpart  for  every  value  of  a in  the  open  interval  (0,  1). 

The  particular  case  a = 0 , which  is  clearly  a limiting  case  of 

non-elliptic  behavior,  requires  separate  comment.  In  this  instance  the 

circle  C is  centered  at  Xj  = 1/2  , x2=  0 and  thus  passes  through  the 

origin.  When  a = 0,  it  turns  out,  the  elliptic  solution  — analogous  to  e 

+ 

of  (3.8)  — can  be  joined  to  the  non-elliptic  solution  corresponding  to  h 
in  (3.  21)  across  the  entire  circle  C in  such  a way  that  both  u and  its 
first  partial  derivatives  remain  continuous.  Hence  the  final  solution  thus 
obtained  is  free  of  shocks.  It  should  be  mentioned  that  the  finite -crack 
problem  (1.  18)  to  (1. 21)  for  a material  obeying  (5.  7)  with  a = 0 is  closely 
related  to  a notch-problem  for  an  elastic-perfectly  plastic  body  in  anti- 
plane shear,  solved  by  Hult  and  McClintock  [7],  who  also  employed  the 
hodograph  method.  In  fact,  for  a notch  of  zero  opening  angle  (crack)  and 
a sufficiently  small  loading  at  infinity,  the  elastic-plastic  interface 
established  in  [7]  — when  referred  to  the  scaling  (2.  6)  — passes  over  into 
the  circle  C that  separates  the  elliptic  from  the  non-elliptic  domain  in 
the  small-scale  nonlinear  crack  problem  corresponding  to  a = 0 in  (5.7). 

An  additional  observation  pertaining  to  the  special  choice  of  cr  = 0 
concerns  the  solution  deduced  in  [2]  to  the  small-scale  crack  problem 
for  the  limiting  elliptic  case  n=  1/2  (see  Figure  2)  of  a power-law 
material,  governed  by  (1.30).  For  such  a material*. 


*See  (6.  52)  in  [2]  and  refer  to  the  scaling  (2.6)  in  the  present  paper. 
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lim  u(r,  9) 

■?-*o 


r 


■u/ (-TTS  6<  -tt/2) 

(i  sin  9 /t^  (- it/2s  0 s tt/2) 

u/ Tco  (tr  / 2 s 0 £ it  ) , 


(5.9) 


where  TQo  = (j/v/Tb  is  the  ultimate  shear  stress.  One  finds  precisely  the 
same  result  for  a simple-shear  response  (5.7)  with  a = 0,  provided 
in  (5.9)  is  taken  to  be  n,  as  is  appropriate.  This  agreement  reflects 
the  fact  that  the  behavior  of  the  solution  at_  the  crack-tip  (as  r-*0)  in 
both  the  limiting  elliptic  case  of  a power -law  material  and  the  limiting 
non-elliptic  case  corresponding  to  a = 0 in  (5.7)  depends  merely  on  the 
ultimate  shear  stress  in  simple  shear. 

An  alternative  modification  of  (1.32)  which  is  of  interest  consists  in 

requiring  that  T(k)  decline  to  a positive  value  t — rather  than  to  zero. 

oo 

as  in  Figure  3.  In  this  event  the  final  solution  again  involves  shocks,  but 

these  are  no  longer  tangential  to  the  Xj-axis  at  the  crack-tip:  instead, 

the  two  shocks  are  now  found  to  subtend  a positive  angle,  which  depends 

on  the  ratio  t /u. 

oo 

Another  study  of  crack  problems  in  which  the  appropriate  differential 
equations  of  equilibrium  are  not  elliptic  is  that  included  in  the  paper  of 
Hutchinson  [8],  who  investigates  the  perfectly  plastic  stress  field  near 
the  tip  of  a crack  in  a thin  sheet  subject  to  tension  at  right  angles  to  the 
crack.  A stress  discontinuity  arises  in  the  results  of  the  analysis  in  [8]. 

In  conclusion  we  remark  that  the  hodograph  scheme  used  in  this 
investigation  is  no  longer  available  in  connection  with  the  physically  more 
important  small-scale  nonlinear  crack  problem  associated  with  a crack 
in  a tension  field  under  conditions  of  plane  strain  and  for  a material 
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capable  of  sustaining  a loss  of  ellipticity.  In  preparation  for  an  attack 
on  this  much  more  difficult  fourth-order  problem,  it  would  be  instructive 
to  find  an  alternative  approach  to  the  crack  problem  treated  here  that  does 
not  rely  on  the  hodograph  method. 
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FIGURE  3.  RESPONSE  CURVE  IN  SIMPLE  SHEAR 
FOR  THE  SPECIAL  PILOT  MATERIAL 
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(c)  u = h (x | , x2)  on#  (HYPERBOLIC  SOLUTION) 

FIGURE  5.  PILOT  MATERIAL:  ELLIPTIC  AND  HYPERBOLIC  SOLUTIONS 


FIGURE  6.  SOLUTION  OF  SMALL-SCALE  NONLINEAR  CRACK  PROBLEM: 
ELLIPTIC  AND  HYPERBOLIC  DOMAINS 


FIGURE  7 ANGULAR  DEPENDENCE  OF  DISPLACEMENT  u 
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FIGURE  8.  ANGULAR  DEPENDENCE  OF  SHEAR  STRESS  T, 


FIGURE  10.  ANGULAR  DEPENDENCE  OF  AXIAL  STRESS  t, 
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